Toward an effective centrality trigger in pp collisions at LHC 
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We investigate the impact of very strong small x gluon fields in colliding nucleons at LHC energies 
on the interaction of valence quarks. We find that in the range of small impact parameters, which 
contribute significantly to the production of heavy new particles, several of the valence quarks 
receive large transverse momenta, exceeding 1 GeV/c. This results in a suppression of leading 
baryon production and consequently in an additional energy flow to smaller rapidities. We suggest 
several triggers for centrality in pp collisions which allow one to study the propagation of partons 
through gluon fields of a strength comparable to the ones encountered in heavy ion collisions at the 
LHC. 

PACS numbers: 12.20.-m, 13.85.Hd 



It is now widely realized that a fast increase of gluon 
densities at small x should result in the onset of the black 
disk regime (BDR) of interactions of partons up to rather 
large virtualities (transverse momenta). In the BDR a 
parton gets a large and with energy increasing trans- 
verse momentum when passing through the gluon field 
(~ Qs m the color glass condensate models, for a recent 
review see In the BDR regime, the parton also loses 
a finite fraction of its energy, which allows one to explain 
the pattern of leading pion production observed at RHIC 
in deuteron-gold collisions [2| . Gluon densities in the nu- 
cleon, at small transverse distances p from the center, 
are comparable to those in nuclei at small p for the same 
energies. In particular for p = 0, and A ~ 200, 
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where po is the nuclear matter density. Here we took 
an exponential shape for the gluon distribution in the 
nucleoli corresponding to the two gluon form factor of 
cxp(— Bt/2) (a dipole fit to the two gluon form factor 
gives an even larger ratio) and neglected the leading twist 
gluon shadowing which suppresses Ga{p) by a factor of 
7a(p). Since snn for pp and for heavy ion collisions at 
the LHC differs by a factor of 6 and the gluon density 
depends on s as s n ,n ~ 0.2 for the virtualities of a few 
GeV 2 , we find that the gluon densities in the two cases 
are similar. Hence up to p ~ 0.6 fm the gluon density 
in a nucleon is comparable to the average gluon nuclear 
density in heavy ion collisions (and larger than the gluon 
densities in heavy ion collisions at RHIC). At larger p the 
nucleon gluon density decreases rapidly. 

There are two challenges: (i) to device a trigger which 
would select pp collisions at small impact parameters 
((b) < (p) ~ 0.6 fm) in order to investigate high gluon 
density effects in pp collisions at densities comparable to 



AA collisions, and (ii) to account for high gluon density 
effects in central pp collisions which give the dominant 
contribution to the production of heavy particles like the 
Higgs, or in SUSY. 

Since the smallest x-i gluons are resolved by large x 
partons in the projectile regions - X2 ~ 4fcf jxs ~ I0~ 7 
for x ~ 0.3 and k t = IGeV/c - we expect the strongest 
effects of the high density gluon fields for the interac- 
tion of valence quarks which should influence the struc- 
ture of the final states in the fragmentation region. Even 
larger effects are present for the interaction of the lead- 
ing gluons due to a factor of 9/4 stronger interaction in 
the gluon-gluon channel. These effects are likely to be 
manifested at smaller rapidities since gluon distributions 
have a much softer x distribution than quarks. (The cur- 
rent MC models usually neglect hard QCD effects in this 
kinematic region focusing on hard QCD dynamics in the 
central region.) 

For a pp collision at a given impact parameter 6, in- 
dividual partons of one nucleon can pass at a range of 
transverse distances p from the second nucleon and hence 
encounter significantly different local gluon densities (see 
FigQJ. In this paper we analyze the effects of the va- 
lence quark interaction with small x gluon fields taking 
into account the geometry of the collisions. This will 
allow us to determine how frequently valence quarks in 
pp collisions at different impact parameters 6, experience 
hard collisions in which they obtain a large transverse 
momentum. Based on this study we propose a series 
of centrality triggers which allow to select collisions at 
much smaller impact parameters than in generic inelas- 
tic events and hence will provide an opportunity to study 
the high gluon field effects in pp collisions. We also sug- 
gest that the pp collisions leading to production of new 
particles like the Higgs boson should be accompanied by a 
significantly stronger flow of energy from the fragmcnta- 
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FIG. 1: Schematic view of the collision geometry. 



tion regions to smaller rapidities than in generic inelastic 
collisions. 

Description of the model. To model the fragmentation 
region in pp collisions we take a simple model for the three 
quark wave function with the distribution of quarks over 
transverse distance from the center given by exp(— Apf ) 
with < p 2 >^ 0.3/m 2 matched to describe the distribu- 
tion of the valence quarks as given by the axial nucleon 
form factor. Accordingly, the event generator produces 
the values of pi for three quarks which are not correlated. 
Note that one does not expect a very strong correlation 
between p's due to the presence of additional partons in 
the wave function (gluons, qq pairs). Nevertheless we 
checked that a requirement | $D i=1 pi\ < 0-1 f m does not 
change results noticeably. Hence we neglect possible cor- 
relations in p between valence quarks. We also assume 
that there are no significant transverse correlations be- 
tween small x (x ~ 10~ 5 ) partons. This assumption is 
based on the presence of diffusion in p in the small x 
evolution which should wash away whatever correlations 
may be present at large x (x > 0.1). (For a discussion of 
the evidence for such correlations see Q.) 

When computing the momentum fractions of the 
quarks, we need to know the virtuality at which the 
quarks are resolved. Since the latter quantity is not 
known beforehand, we generate x B ,i and pi from dx/x — 
const, and dp = const, distributions. The selection ac- 
cording to the structure functions and the form factor is 
done in the end, after specifying Q 2 , via rejection. For 
a given impact parameter b and relative positions of the 
quarks pi in the projectile, we estimate the density within 
the color glass condensate approach. So we need to know 
the saturation scale in the target for the three valence 
quarks. 



Q 2 = Q 2 s {x A ,\b- 
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with x a = Q 2 / (sxb). The saturation momentum is pa- 
rameterized as 



\{x A ,p) = Qs, ( ) F g (x A ,p;Q 2 )/c F , (>!) 
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FIG. 2: Probability for the different classes of events with n 
quarks struck at a given impact parameter b. 



where cf normalizes the density. We choose xq = 0.01, 
Q 2 = 0.6 GeV 2 and cf = F g (xo,0;Q 2 ) such that 
the saturation momentum in the center of the target at 
x A = xq is just Q 2 . The implicit definition for the sat- 
uration scale in eq. ([3]) is solved by a simple iteration, 
the expression converges after a few steps. Finally, the 
whole configuration is accepted with the probability 



p ~ pF g {x Bl p;Q 2 s )x B fGYiy{xB,Q 2 s ) 



(4) 



where x/grv are standard GRV structure functions of 
the proton, and the two-gluon form factor at high mo- 
mentum fraction xb describes the spatial distribution of 
the valence quarks. The actual transverse momentum 
kick is then drawn from the distribution 0, [B| 
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We conservatively considered only the case when the 
BDR is reached for Q s > 1 GcV/c and counted only 
quark interactions in which the quark received a trans- 
verse momentum k t > 0.75GeV/c. The reason for such a 
cut is that for such momenta, the probability to form a 
nucleon with large longitudinal momentum is suppressed, 
as a minimum, by the square of the nucleon form factor 
Fjy(fct). In the BDR a quark not only gets a large trans- 
verse momentum but also loses a finite fraction of its mo- 
mentum - we neglected this effect in our analysis since it 
is likely to be on the scale of a 10% energy loss and would 
not affect our qualitative discussion. We also assumed in 
the analysis that the gluon densities at small x are not 
fluctuating significantly for a impact parameters where 
the BDR holds (this is an assumption usually made in 
the saturation models where the the small gluon field is 
treated as nearly classical field). Neglecting fluctuations 
seems safe for the small x in our analysis, though for 
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FIG. 3: dn/dxF distribution of baryons. FIG. 4: Probability that no quarks were struck (solid curve 

as compared to the probability of absence of generic inelastic 
interaction (dashed curve) as a function of b. 



larger x relevant at the Tevatron collider they may be 
more important. 

As a first result, the probabilities Pi for i quarks to in- 
teract in the BDR regime (number of "wounded" quarks, 
N w ) are presented in Fig. [2] as a function of the impact 
parameter b. Note that for small b interactions N w > 2 
dominate. In such interactions each quark fragments in- 
dependently, so the spectra for N w — 2 and N w = 3 
should be rather similar and shifted to much smaller xf 
than in soft interactions where the spectra of nucleons are 
known to be flat in xf in a wide range of xf ■ At the same 
time for N w = 1, the dominant process is likely to be 
diquark fragmentation which is similar to the soft mech- 
anism. Hence one expects two distinctive shapes of xp 
distributions - one with leading baryons (N w = 0,1) and 
another with suppressed forward scattering (N w = 2, 3). 
We implemented the fragmentation of the system pro- 
duced in the first stage by constructing strings which 
decay using the LUND method. There are always two 
strings, drawn between a quark and a diquark from the 
interacting particles. When a quark of the diquark re- 
ceives a high transverse momentum, the diquark becomes 
a system of two quarks and a junction. This has the nice 
property that one recovers the diquark when the invari- 
ant mass between the two quarks is small. The results 
are plotted in Fig. [3] and are in good agreement with the 
qualitative expectations discussed above. 

The total probability of an inelastic interaction at 
given impact parameter can be expressed through the 
impact factor T(b), which is the Fourier transform of 
the elastic amplitude, as 1 — |1 — r(6)| 2 . On the other 
hand in our model we calculate the probability of inelas- 
tic interactions due to the hard interaction of i quarks 
of one of the nucleons, Pi. Since in the discussed ap- 
proximation there is no correlation between the interac- 
tions in two fragmentation regions for a fixed b the prob- 
ability that none of six valence quarks will get a kick is 



1 — (1 — Pi — Pi — P3) 2 . For a consistency of the model 
we need that 

i-|i-r(6)| 2 >i-(i- p- P2-P3) 2 , (6) 

or 

r(6) >Pi+P 2 + P 3 . (7) 

Wc find that this condition is satisfied for all b (within 
the uncertainties of the current models for pp elastic scat- 
tering at LHC) - see Fig. 0] indicating that hard collisions 
lead to saturation of the black limit for V for small b at 
LHC energies, while at large b hard interactions give a 
very small contribution. Hence the model provides a dy- 
namic explanation of T(b ~ 0) = 1 and leads to expansion 
of the r(6) = 1 region with increase of energy. 

Very strong dependence of N w on b and a strong cor- 
relation of N w with the multiplicity of leading baryons 
allows one to determine the effectiveness of a central- 
ity trigger based on a veto for the production of lead- 
ing baryons with x > Xt r as a function of Xt r - We find 
than an optimal value of Xt r is ~ 0.1. Current config- 
urations of several LHC detectors allow to veto neutron 
production in this x-range. TOTEM, in addition, al- 
lows to veto production of protons with xf > 0.8. Since 
neutron and proton multiplicities are similar, a one side 
veto for production of both charged and neutral baryons 
leads approximately to the same result as a two side veto 
for neutron production. Accordingly we will give results 
both for single side veto and for two side veto for both 
neutral and charged baryons (understanding that the full 
implementation of the latter option would require certain 
upgrades of the detectors some of which are currently un- 
der discussion). 

The results of the calculations are presented in Fig. [5] 
together with the distribution over b for generic inelastic 
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FIG. 5: Impact parameter distributions for inelastic events, 
the dijet trigger and single and double sided veto-trigger (no 
baryon in the region xf > 0.1). 
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FIG. 6: The combination of dijet and veto trigger gives the 
best constraints on central events in pp-collisions. 



events and the central dijet trigger [6|. We see that the 
single side veto trigger leads to a centrality similar to that 
of a the dijet trigger, while a double side veto leads to the 
most narrow distribution in b. An easy way to check this 
expectation would be to compare other characteristics of 
events due to these triggers - one expects for example 
a progressive increase of the central multiplicity with a 
decrease of average b. 

The most narrow distributions can be achieved by 
combining two types of the triggers - veto and dijet - 
Fig. O We find that in this case we reach the limit 

1 /2 

that { Ptr ) = ((p 2 ) + (b 2 )) ' beco mes comparable to (p) 
which is the smallest possible average (p) for pp or DIS 
collisions. 

In line with discussion of Eq. Q] we expect that for our 
best centrality triggers average gluon densities will be 
comparable to those in the heavy ion collisions. Indeed, 



we estimate the average gluon density encountered by 
the leading partons, and find it larger than the one en- 
countered at central impact parameters in pA collisions 
at RHIC where BDR holds at least for p t < 1.5GeV/c. 
Further, it is about \ of the average gluon density in 
heavy ion collisions at LHC. (However dispersion in the 
value of Q s is much larger in the pp case). 



We also estimated the effect of a dijet trigger (which is 
the same selection of impact parameters as for production 
of heavy new particles). We find that in 16% of these 
events at least in one fragmentation region no leading 
baryons would be present with xp > 0.1 as compared to a 
generic case where this fraction is 8%. This should result 
in a significantly larger flow of energy to central rapidities 
and strong fluctuations of this flow on an event by event 
basis. Since the probability of interaction for gluons with 
x > 10~ 2 is at least as large as for quarks with x > 0.1 we 
expect that many of the gluons with these x will also lose 
coherence with the low transverse momentum partons 
and will fragment independently further contributing to 
the energy flow to central rapidities. It is important to 
investigate how these effect would impact on the various 
triggers suggested for the searches of the new particles. 

In conclusion, we have demonstrated that it will be 
feasible to trigger on very central pp collisions at the 
LHC. The high gluon densities encountered by the col- 
liding partons in such collisions facilitate studies of novel 
QCD effects complementary to heavy ion collisions, such 
as the p t broadening of the Z-boson distribution as com- 
pared to generic events, multiplicity fluctuations, elliptic 
flow, etc. 
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